Turbulent Disks are Never Stable: Fragmentation and Turbulence-Promoted Planet Formation

A fundamental assumption in our understanding of disks is that when the Toomre Q>>1, the disk is stable against fragmentation into self-gravitating objects (and so cannot form planets via direct collapse). But if disks are turbulent, this neglects a spectrum of stochastic density fluctuations that can produce rare, high-density mass concentrations. Here, we use a recently-developed analytic framework to predict the statistics of these fluctuations, i.e. the rate of fragmentation and mass spectrum of fragments formed in a turbulent Keplerian disk. Turbulent disks are never completely stable: we calculate the (always finite) probability of forming self-gravitating structures via stochastic turbulent density fluctuations in such disks. Modest sub-sonic turbulence above Mach number ~0.1 can produce a few stochastic fragmentation or 'direct collapse' events over ~Myr timescales, even if Q>>1 and cooling is slow (t_cool>>t_orbit). In trans-sonic turbulence this extends to Q~100. We derive the true Q-criterion needed to suppress such events, which scales exponentially with Mach number. We specify to turbulence driven by MRI, convection, or spiral waves, and derive equivalent criteria in terms of Q and the cooling time. Cooling times >~50*t_dyn may be required to completely suppress fragmentation. These gravoturbulent events produce mass spectra peaked near ~M_disk*(Q*M_disk/M_star)^2 (rocky-to-giant planet masses, increasing with distance from the star). We apply this to protoplanetary disk models and show that even minimum mass solar nebulae could experience stochastic collapse events, provided a source of turbulence.Comment: 15 pages, 5 figures (+appendix), accepted to ApJ (added clarifications and discussion to match accepted version

Bridging the Gap: Observations and Theory of Star Formation Meet on Large and Small Scales

The drive to understand galaxy formation and evolution over the lifetime of the universe has justified vast space-based and ground-based telescope facilities, as well as the development of new technologies. The details of star formation, and the modes by which that activity couples to the broader galactic environment, occurs on small spatial scales. These scales can only be traced with great sophistication in the local universe, as witnessed by observations using the Spitzer Space Telescope, the Herschel Space Observatory, the Stratospheric Observatory for Infrared Astronomy (SOFIA), and the Atacama Large Millimeter Array (ALMA). A critical realization of the last decade, however, is that the large scales and small scales are strongly coupled, and cannot be treated in isolation. At the same time, the scales studied in the local universe are "sub-grid" for the purpose of cosmological simulations that make valiant efforts to include the physics of star formation and its feedback to the local environment. The fundamental uncertainties in how this sub-grid physics is incorporated into the larger picture are by far the greatest limitation in understanding galaxy formation and evolution. The main goal of cosmological simulations is to understand the formation and evolution of the universe over a wide range of scales going from the Hubble volume to sub-light-year scales within galaxies. However, due to computational limitations, the physics of star formation in galaxies and the effects this process has on the formation and evolution of galaxies are often greatly simplified. For instance, accounting for the effect of stellar feedback in cosmological simulations is a crucial factor in galaxy formation and evolution. Without stellar feedback in galaxies, the gas would rapidly cool and collapse, converting all available gas into stars within a dynamical time. This consequence is in sharp conflict with observations—there are vastly fewer stars in our universe than the models would predict. Until recently, numerical simulations have been unable to regulate star formation efficiently enough to reproduce observations, with many consequences: models with too many stars also predict they form at the wrong times in the universe’s history, that there are the wrong abundances of heavy elements, that galaxies look nothing like the Milky Way, and even that the observable properties of dark matter are (apparently) discordant with new precision-cosmology measurements. All of these problems may, in fact, be due to the fundamental problem of understanding how small and large scales interact

Galaxy disks do not need to survive in the L-CDM paradigm: the galaxy merger rate out to z~1.5 from morpho-kinematic data

About two-thirds of present-day, large galaxies are spirals such as the Milky Way or Andromeda, but the way their thin rotating disks formed remains uncertain. Observations have revealed that half of their progenitors, six billion years ago, had peculiar morphologies and/or kinematics, which exclude them from the Hubble sequence. Major mergers, i.e., fusions between galaxies of similar mass, are found to be the likeliest driver for such strong peculiarities. However, thin disks are fragile and easily destroyed by such violent collisions, which creates a critical tension between the observed fraction of thin disks and their survival within the L-CDM paradigm. Here we show that the observed high occurrence of mergers amongst their progenitors is only apparent and is resolved when using morpho-kinematic observations which are sensitive to all the phases of the merging process. This provides an original way of narrowing down observational estimates of the galaxy merger rate and leads to a perfect match with predictions by state-of-the-art L-CDM semi-empirical models with no particular fine-tuning needed. These results imply that half of local thin disks do not survive but are actually rebuilt after a gas-rich major merger occurring in the past nine billion years, i.e., two-thirds of the lifetime of the Universe. This emphasizes the need to study how thin disks can form in halos with a more active merger history than previously considered, and to investigate what is the origin of the gas reservoir from which local disks would reform.Comment: 19 pages, 7 figures, 2 tables. Accepted in ApJ. V2 to match proof corrections and added reference

On The Nature of Variations in the Measured Star Formation Efficiency of Molecular Clouds

Measurements of the star formation efficiency (SFE) of giant molecular clouds (GMCs) in the Milky Way generally show a large scatter, which could be intrinsic or observational. We use magnetohydrodynamic simulations of GMCs (including feedback) to forward-model the relationship between the true GMC SFE and observational proxies. We show that individual GMCs trace broad ranges of observed SFE throughout collapse, star formation, and disruption. Low measured SFEs (<<1%) are "real" but correspond to early stages, the true "per-freefall" SFE where most stars actually form can be much larger. Very high (>>10%) values are often artificially enhanced by rapid gas dispersal. Simulations including stellar feedback reproduce observed GMC-scale SFEs, but simulations without feedback produce 20x larger SFEs. Radiative feedback dominates among mechanisms simulated. An anticorrelation of SFE with cloud mass is shown to be an observational artifact. We also explore individual dense "clumps" within GMCs and show that (with feedback) their bulk properties agree well with observations. Predicted SFEs within the dense clumps are ~2x larger than observed, possibly indicating physics other than feedback from massive (main sequence) stars is needed to regulate their collapse.Comment: Fixed typo in the arXiv abstrac

Microjansky sources at 1.4 GHz

We present a deep 1.4 GHz survey made with the Australia Telescope Compact Array (ATCA), having a background RMS of 9 microJy near the image phase centre, up to 25 microJy at the edge of a 50' field of view. Over 770 radio sources brighter than 45 microJy have been catalogued in the field. The differential source counts in the deep field provide tentative support for the growing evidence that the microjansky radio population exhibits significantly higher clustering than found at higher flux density cutoffs. The optical identification rate on CCD images is approximately 50% to R=22.5, and the optical counterparts of the faintest radio sources appear to be mainly single galaxies close to this optical magnitude limit.Comment: 6 pages, 4 figures, accepted by ApJ Letters 4 May 199

The Origin and Evolution of the Galaxy Mass-Metallicity Relation

We use high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environment (FIRE) project to study the galaxy mass-metallicity relations (MZR) from z=0-6. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback. The simulations cover halo masses Mhalo=10^9-10^13 Msun and stellar mass Mstar=10^4-10^11 Msun at z=0 and have been shown to produce many observed galaxy properties from z=0-6. For the first time, our simulations agree reasonably well with the observed mass-metallicity relations at z=0-3 for a broad range of galaxy masses. We predict the evolution of the MZR from z=0-6 as log(Zgas/Zsun)=12+log(O/H)-9.0=0.35[log(Mstar/Msun)-10]+0.93 exp(-0.43 z)-1.05 and log(Zstar/Zsun)=[Fe/H]-0.2=0.40[log(Mstar/Msun)-10]+0.67 exp(-0.50 z)-1.04, for gas-phase and stellar metallicity, respectively. Our simulations suggest that the evolution of MZR is associated with the evolution of stellar/gas mass fractions at different redshifts, indicating the existence of a universal metallicity relation between stellar mass, gas mass, and metallicities. In our simulations, galaxies above Mstar=10^6 Msun are able to retain a large fraction of their metals inside the halo, because metal-rich winds fail to escape completely and are recycled into the galaxy. This resolves a long-standing discrepancy between "sub-grid" wind models (and semi-analytic models) and observations, where common sub-grid models cannot simultaneously reproduce the MZR and the stellar mass functions.Comment: 17 pages, 14 figures, re-submitted to MNRAS after revisions on referee comment